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Abstract. One of the topical problems of contemporary physics is a possible variability 
of the fundamental constants. Here we consider possible variability of two dimensionless 
constants which are most important for calculation of atomic and molecular spectra (in 
particular, the X-ray ones): the fine-structure constant a = e 2 /tic and the proton-to- 
electron mass ratio // = m p /m e . Values of the physical constants in the early epochs 
are estimated directly from observations of quasars - the most powerful sources of 
radiation, whose spectra were formed when the Universe was several times younger 
than now. A critical analysis of the available results leads to the conclusion that 
present-day data do not reveal any statistically significant evidence for variations of 
the fundamental constants under study. The most reliable upper limits to possible 
variation rates at the 95% confidence level, obtained in our work, read: 

\a/a\ < 1.4 x 1(T 14 yr" 1 , \fi/fi\ < 1.5 x 1(T 14 yr" 1 

on the average over the last 10 10 yr. 



INTRODUCTION 

Contemporary theories (SUSY GUT, superstring and others) not only predict 
the dependence of fundamental physical constants on energy 1 , but also have cos- 
mological solutions in which low-energy values of these constants vary with the 
cosmological time. The predicted variation at the present epoch is small but non- 
zero, and it depends on theoretical model. In particular, Damour and Polyakov 
[1] have developed a modern version of the string theory, whose parameters could 
be determined from cosmological variations of the coupling constants and hadron- 
to-electron mass ratios. Clearly, a discovery of these variations would be a great 
step in our understanding of Nature. Even a reliable upper bound on a possible 
variation rate of a fundamental constant presents a valuable tool for selecting viable 
theoretical models. 



^ The prediction of the theory that the fundamental constants depend on the energy of interac- 
tion has been confirmed in experiment. In this paper, we consider only the space-time variability 
of their low-energy limits. 



Historically, a hypothesis that the fundamental constants may depend on the 
cosmological time t (that is the age of the Universe) was first discussed by Milne [2] 
and Dirac [3]. The latter author proposed his famous "large- number hypothesis" 
and suggested that the gravitational constant was directly proportional to t. Later 
the variability of fundamental constants was analyzed, using different arguments, 
by Gamow [4] , Dyson [5] , and others. The interest in the problem has been revived 
due to recent major achievements in GUT and Superstring models (e.g., [1]). 

Presently, the fundamental constants are being measured with a relative error 
of ~ 10 -8 . These measurements obviously rule out considerable variations of the 
constants on a short time scale, but do not exclude their changes over the lifetime 
of the Universe, ~ 1.5 x 10 10 years. Moreover, one cannot rule out the possibility 
that the constants differ in widely separated regions of the Universe; this could be 
disproved only by astrophysical observations and different kinds of experiments. 

Laboratory experiments cannot trace possible variation of a fundamental con- 
stant during the entire history of the Universe. Fortunately, Nature has provided 
us with a tool for direct measuring the physical constants in the early epochs. This 
tool is based on observations of quasars, the most powerful sources of radiation. 
Many quasars belong to most distant objects we can observe. Light from the dis- 
tant quasars travels to us about 10 10 years. This means that the quasar spectra 
registered now were formed ~ 10 10 years ago. The wavelengths of the lines ob- 
served in these spectra (A Q b s ) increase compared to their laboratory values (Ai a b) 
in proportion A Q b s = Ai a b(l + z), where the cosmological redshift z can be used to 
determine the age of the Universe at the line-formation epoch. In some cases, the 
redshift is as high as z ~ 3—5, so that the intrinsically far-ultraviolet lines are reg- 
istered in the visible range. The examples are demonstrated in Fig. 1. Analysing 
these spectra we may study the epoch when the Universe was several times younger 
than now. 

Here we review briefly the studies of the space-time variability of the fine- 
structure constant a and the proton-to-electron mass ratio ji. 



FINE-STRUCTURE CONSTANT 

Various tests of the fundamental constant variability differ in space-time regions 
of the Universe which they cover. Local tests relate to the values of constants on 
the Earth and in the Solar system. In particular, laboratory tests infer the possible 
variation of certain combinations of constants "here and now" from comparison of 
different frequency standards. Geophysical tests impose constraints on combina- 
tions of fundamental constants over the past history of the Solar system, although 
most of these constraints are very indirect. In contrast, astrophysical tests allows 
one to "measure" the values of fundamental constants in distant areas of the early 
Universe. 
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FIGURE 1. Portions of quasar spectra which show the absorption lines of H2 and Siiv with 
large rcdshifts, z — (A Q b s — Ai a b)/Ai a b ~ 2.8. The lower horizontal axis gives the wavelengths in 
the observer's frame (A bs) and the upper axis gives the wavelengths in the quasar's frame Ai a b (in 
A), (a) The spectrum (thin line) of the quasar PKS 0528-250, obtained with the 4- meter CTIO 
telescope (Chile), containing H 2 lines which belong to the L 4-0 branch of the spectrum; thick 
line plots the spectral fit. (b) The spectrum of the quasar HS 1946+76, obtained with the 6-meter 
SAO telescope (Russia), containing Siiv doublet lines which correspond to the 2 S!/ 2 ^ 2 P3/2 and 
2 Si/2 ^ 2 Pi/2 transitions. 



Local tests 



Laboratory experiments 

There were a number of laboratory experimets aimed at detection of trends of 
the fundamental constants with time by comparison of frequency standards which 
have different dependences on the constants. We mention only two of the published 
experiments. 

Comparison of H-masers with Cs-clocks during 427 days revealed a relative (H- 
Cs) frequency drift with a rate 1.5 x 10~ 16 per day, while the rates of (H-H) and 
(Cs-Cs) drifts (i.e., the drifts between identical standards, used to control their 
stability) were less than 1 x 10~ 16 per day [6]. A similar result was found in 
comparison of a Hg + -clock with a H-maser during 140 days [7]: the rate of the 
relative frequency drift was less than (2 ± 1) x 10~ 16 per day. 

Such a drift is treated as a consequence of a difference in the long-term stability 



of different atomic clocks. In principle, however, it may be caused by variation of a. 
That is why it gives an upper limit to the a variation [7]: \a/a\ < 3.7 x ICT 14 yr _1 . 

Geophysical tests 

The strongest bound on the possible time-variation rate of a was derived in 
1976 by Shlyakhter [8], and recently, from a more detailed analysis, by Damour 
and Dyson [9], who obtained \a/a\ < 0.7 x ICT 16 yr _1 , The analysis was based 
on measurements of isotope ratios in the Oklo site in Gabon, where a unique 
natural uranium nuclear fission reactor had operated 1.8 billion years ago. The 
isotope ratios of samarium produced in this reactor by the neutron capture reaction 
149 Sm+n ^ 150 Sm+7 would be completely different, if the energy of the nuclear 
resonance responsible for this capture were shifted at least by 0.1 eV. 

Another strong bound, \a/a\ < 5 x 10 -15 yr" 1 , was obtained by Dyson [5] from 
an isotopic analysis of natural radioactive decay products in meteorites. 

A weak point of these tests is their dependence on the model of the phenomenon, 
fairly complex, involving many physical effects. For instance, Damour and Dyson 
[9] estimated possible shift of the above-mentioned resonance due to the a variation, 
assuming that the Coulomb energy of the excited state of 150 Sm*, responsible for 
the resonance, is not less than the Coulomb energy of the ground state of 150 Sm. 
In absence of experimental data on the nuclear state in question, this assumption 
is not justified, since heavy excited nuclei often have Coulomb energies smaller 
than those for their ground states [10]. Furthermore, a correlation between the 
constants of strong and electroweak interactions (which is likely in the frame of 
modern theory) might lead to further softening of the mentioned bounds by 100- 
fold, to \a/a\ < 5 x 10~ 15 yr _1 , as noted by Sisterna and Vucetich [11]. 

In addition, the local tests cannot be extended to distant space regions and to 
the early Universe, since the law of possible space-time variation of a is unknown 
a priory. It is the extragalactic astronomy that allows us to study these remote 
regions of spacetime, in particular the regions which were causally disconnected at 
the epoch of formation of the observed absorption spectra. 

Astrophysical tests 

To find out whether a changed over the cosmological time, we have studied the 
fine splitting of the doublet lines of Si IV, C IV, Mg II and other ions, observed in 
the spectra of distant quasars. According to quantum electrodynamics, the relative 
splitting of these lines 5 A/ A is proportional to a 2 (neglecting very small corrections). 
Consequently, if a changed with time, then 5\/\ would depend on the cosmological 
redshift z. This method of measuring a in distant regions of the Universe had been 
first suggested by Savedoff [12] and was used later by other authors. For instance, 
Wolfe et al. [13] derived an estimate \a/a\ < 4 x 10~ 12 yr _1 from an observation of 
the Mgll absorption doublet at z = 0.524. 



TABLE 1. Variation of a value estimated from 
rcdshiftcd Si iv fine- splitting doublets. 
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An approximate formula which relates a deviation of a at redshift z from its cur- 
rent value, Aa z , with measured S\/\ in the extragalactic spectra and in laboratory 
reads 



Aa z f« — 
2 



(<WA) 2 
(5A/A) 



where c r ~ 1 takes into account radiation corrections [14]: for instance, for Siiv 
c r w 0.9. 

Many high-quality quasar spectra measured in the last decade have enabled us to 
significantly increase the accuracy of determination of <5A/A at large z. An example 
of the spectra observed is shown in Fig. I. For the present report, we have selected 
the results of high-resolution observations [15-17], most suitable for an analysis of 
a variation. The values of Aa/a calculated from these data according to Eq. (1) 
are given in Table 1. 

As a result, we obtain a new estimate of the possible deviation of the fine- 
structure constant at z — 2-4 from its present [z = 0) value: 

Aa/a = (-4.6 ± 4.3 [stat] ± 1.4 [syst]) x 1(T 5 , (2) 

where the statistical error is obtained from the scatter of astronomical data (at large 
z) and the systematic one is estimated from the uncertainty of the fine splitting 
measurement in the laboratory [18,19] (at z — 0, which serves as the reference point 
for the estimation of Aa). The corresponding upper limit of the a variation rate 
averaged over ~ 10 10 yr is 



\a/a\ < 1.4 x 1(T 14 yr" 1 (3) 

(at the 95% confidence level). This constraint is much more stringent than those 
obtained from all but one previous astronomical observations. The notable ex- 
ception is presented by Webb et al. [20], who have analysed spectroscopic data of 
similar quality, but estimated a from comparison of Fell and Mgll fine-splitted 
walelengths in extragalactic spectra and in the laboratory. Their result indicates a 
tentative time-variation of a: Aa/a = (—1.9 ± 0.5) x 10~ 5 at z — 1.0-1.6. Note, 
however, two important sources of a possible systematic error which could mimic 
the effect: (a) Fe,ll and Mgll lines used are situated in different orders of the 
echelle-spectra, so relative shifts in calibration of the different orders can simu- 
late the effect of a-variation, and (b) were the relative abundances of Mg isotopes 
changing during the cosmological evolution, the Mg II lines would be subjected to 
an additional ^-dependent shift relative to the Fe II lines, quite sufficient to sim- 
ulate the variation of a (this shift can be easily estimated from recent laboratory 
measurements [21]). In contrast, the method based on the fine splitting of a line 
of the same ion species (Si IV in the above example) is not affected by these two 
uncertainty sources. Thus we believe that the restriction (3) is the most reliable at 
present for the long-term history of the Universe. 

According to our analysis, some theoretical models are inconsistent with observa- 
tions. For example, power laws a oc t n with n — 1, —1/4, and —4/3, published by 
various authors in 1980s, are excluded. Moreover, the Teller-Dyson's hypothesis on 
the logarithmic dependence of a on t [22,5] has also been shown to be inconsistent 
with observations. 

Many regions of formation of the spectral lines, observed at large redshifts in 
different directions in the sky, had been causally disconnected at the epochs of line 
formation. Thus, no information could have been exchanged between these regions 
of the Universe and, in principle, the fundamental constants could be different 
there. However, a separate analysis [23] has shown that a value is the same in 
different directions in the sky within the 3a relative error \Aa/a\ < 3 x 10 -4 . 

PROTON-TO-ELECTRON MASS RATIO 

The dimensionless constant /i = m p /m e approximately equals the ratio of the 
constants of strong interaction g 2 / (he) ~ 14 and electromagnetic interaction a ~ 
1/137.036, where g is the effective coupling constant calculated from the amplitude 
of nucleon-7r-meson scattering at low energy. 

In order to check the cosmological variability of \x we have used high-redshift 
absorption lines of molecular hydrogen H2 in the spectrum of the quasar PKS 
0528-250. This is the first (and, in a sense, unique) high-redshift system of H2 
absorption lines discovered in 1985 [24]. A study of these objects yields information 
of paramount importance on the physical conditions ~ 10 10 years ago. 

A possibility of distinguishing between the cosmological redshift of spectral wave- 
lengths and shifts due to a variation of ji arises from the fact that the electronic, 



TABLE 2. Comparison of wavelengths of clcctron-vibro-rotational 
lines for H 2 , D 2 , and T 2 . 



i MH 2 ) A t (D 2 ) A,(T 2 ) Kj 

L 0-0 R(l) 1108.633 1103.351 1101.021 -8.18 x 10~ 3 

L 0-2 R(l) 1077.697 1081.153 1082.760 +5.35 x 10~ 3 

L 0-9 R(l) 992.013 1015.610 1027.218 +3.80 x 10~ 2 



vibrational, and rotational energies of H 2 each undergo a different dependence on 
the reduced mass of the molecule. Hence comparing ratios of wavelengths Aj of 
various H 2 electron-vibration-rotational lines in a quasar spectrum at some red- 
shift z and in laboratory (at z = 0), we can trace variation of /x. The method had 
been used previously by Foltz et al. [25], whose analysis was corrected later in our 
papers [26,23,27]. In the latter papers, we calculated the sensitivity coefficients Ki 
of the wavelengths \ with respect to possible variation of /i and applied a linear 
regression analysis to the measured redshifts of individual lines Z{ as function of 
An illustration of the wavelength dependences on the mass of the nucleus is 
given in Table 2, where a few resonance wavelengths of hydrogen, deuterium, and 
tritium molecules are listed. One can see that, as the nuclear mass increases, dif- 
ferent wavelengths shift in different directions. More complete tables, as well as 
two algorithms of Ki calculation, are given in Refs. [23,28]. 

Thus, if the proton mass in the epoch of line formation were different from the 
present value, the measured Zi and Ki values would correlate: 

We have performed a z-to-K regression analysis using a modern high-resolution 
spectrum of PKS 0528—250 [28]. Several tens of the H 2 lines have been identified; 
a portion of the spectrum which reveales some of the lines is shown in Fig. 1. The 
redshift estimates for individual absorption lines with their individual errorbars 
are plotted in Fig. 2 against their sensitivity coefficients. The resulting parameter 
estimate and la uncertainty is 

A/V/t = (-11-5 ± 7.6 [stat] ± 1.9 [syst]) x 10" 5 . (5) 

The 2a confidence interval to A/1//1 is 

|A/VH < 2.0 x 10~ 4 . (6) 

Assuming that the age of the Universe is ~ 15 Gyr the redshift of the H 2 absorption 
system z = 2.81080 corresponds to the elapsed time of 13 Gyr (in the standard 
cosmological model). Therefore we arrive at the restriction 

\(i/n\ < 1.5 x 10~ 14 yr- 1 (7) 

on the variation rate of /x, averaged over 90% of the lifetime of the Universe. 
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FIGURE 2. Rcdshift values inferred from an analysis of separate spectral features in an H2 
absorption system in the spectrum of the quasar PKS 0528—250, plotted vs. Ai (/i)-sensitivity 
coefficients Ki. The slanted solid line shows the most probable regression and the dashed ones 
corespond to ±la deviations of the slope. 

CONCLUSIONS 

Despite the theoretical prediction of the time-dependences of fundamental con- 
stants, a statistically significant variation of any of the constants have not been 
reliably detected up to date, according to our point of view substantiated above. 
The upper limits obtained indicate that the constants of electroweak and strong 
interactions did not significantly change over the last 90% of the history of the 
Universe. This shows that more precise measurements and observations and their 
accurate statistical analyses are required in order to detect the expected variations 
of the fundamental constants. 
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